Technical 

Research 

Report 


SYSTEMS 

RESEARCH 

CENTER 


Supported  by  the 
National  Science  Foundation 
Engineering  Research  Center 
Program  (NSFD  CD  8803012), 
Industry  and  the  University 

Range  of  the  /c-Dimensional  Radon  Transform 
in  Real  Hyperbolic  Spaces 


by  C.A.  Berenstein  and  E.C.  Tarabusi 


TR  91-101 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

1991 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1991  to  00-00-1991 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Range  of  the  k-Dimensional  Radon  Transform  in  Real  Hyperbolic  Spaces 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Maryland, Systems  Research  Center, College 

Park, MD, 20742 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

see  report 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF  PAGES 

19 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Range  of  the  fc-Dimensional  Radon  Transform 
in  Real  Hyperbolic  Spaces 

Carlos  A.  BERENSTEIN*  and  Enrico  CASADIO  TARABUSH 


Abstract.  Characterizations  of  the  range  of  the  totally  geodesic  k-dimensional  Radon  trans¬ 
form  on  the  n-dimensional  hyperbolic  space  are  given  both  in  terms  of  moment  conditions 
and  as  the  kernel  of  a  differential  operator. 


1.  Introduction 

The  recent  interest  on  the  ^-dimensional  totally  geodesic  Radon  transform  R  in  the  real 
7t- dimensional  hyperbolic  space  Hn  has  arisen  in  part  due  to  its  applications  to  Electrical 
Impedance  Tomography  (EIT).  This  transform  was  introduced  by  Helgason  in  [Ell],  who 
found  two  kinds  of  inversion  formulas,  the  former  valid  for  k  even  [Hi],  [H2],  the  latter  for 
any  k  [H3].  Since  the  approximate  inversion  algorithm  for  EIT  proposed  by  Barber  and 
Brown  [BBl],  [BB2],  [SV],  amounts  to  using  the  backprojection  R*  as  an  approximate 
inverse,  the  authors  searched  for  a  filtered  backprojection  inversion  [BC].  The  crucial  step 
is  Helgason’s  observation  that  R* R  acts  on  the  space  <S(Hn)  as  a  convolution  operator 
with  a  radial  function,  namely,  up  to  a  multiplicative  constant,  sinh k~n  r.  Using  symbolic 
calculus  one  finds  the  algorithm  proposed  in  [BC],  which  has  the  form 

Pn,k(&H)Sn,k  *  R*R  =  I, 

where  pUjk  is  a  polynomial  in  the  Laplace-Beltrami  operator  A h  and  Sn,k  is  a  radial 
integrable  function.  By  analogy  with  the  Euclidean  Radon  transform,  one  expects  that 
this  algorithm  will  perform  better  than  [BBl],  [BB2],  [IC],  [AS],  in  solving  the  EIT 
defining  equations. 

Another  question  that  arises  frequently  in  applications  is  the  characterization  of  the 
ranges  of  the  Radon  transforms.  This  question,  whose  answer  was  long  known  in  the  space 
H(R")  for  the  Euclidean  Radon  transform  [H2,  Corollary  1.2.28]  has  only  recently  been 
completely  settled  in  <S(Rn).  We  have  found  intertwining  operators  between  Euclidean 
and  hyperbolic  Radon  transforms,  hence  the  characterization  of  the  range  of  the  Radon 
transform  in  P(Hn)  and  <S(H")  can  be  settled  very  easily.  Another  consequence  of  the 
existence  of  intertwining  operators  is  that  one  finds  other  inversion  formulas  for  the  hy¬ 
perbolic  Radon  transform.  This  was  inspired  by  the  works  of  Quinto  [Q]  and  Kurusa 
[Ku2]. 

The  authors  would  like  to  thank  W.  M.  Goldman,  E.  L.  Grinberg,  E.  T.  Quinto,  and 
R.  S.  Strichartz  for  informative  conversations. 


^Research  partially  supported  by  grants  NSFDCDR-8803012,  NSFDMS  90-00619,  and  AFOSR  90-0105. 
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2.  Preliminaries 


We  recall  here  the  definitions  and  notation  for  the  Radon  and  Riesz  transforms  in  Hn 
as  given  in  [H2].  We  shall  use  the  ‘conformal  disk’  model  for  Hn,  viz.,  the  open  unit  ball 
Bn  of  Rn  with  the  metric 


dsJ 


4  dx2 


4  £"=1*1 
(i-M2)2  (i-E£,4)! 


where  ||  •  ||  denotes  the  Euclidean  norm  in  Rn.  Such  metric  is  in  fact  conformal  to  the 
Euclidean  one  dx 2  and  has  constant  curvature  —1  (in  some  chapters  of  [H2]  the  curvature 
is  —4).  The  induced  distance  between  x,y  G  Hn  is 


d(x,  y)  =  2  arcsinh 


-  yii 


conversely,  the  Euclidean  norm  of  x  can  be  recovered  by 


11*11  =  tanh  . 

The  geodesics  and  the  totally  geodesic  hypersurfaces  of  Hn  are  arcs  of  circle,  respectively 
spherical  caps,  which  intersect  Sn_1  perpendicularly.  The  spheres  which  are  tangent  to 
S"-1  and  are  contained  in  the  unit  ball  are  called  horocycdes. 

In  geodesic  polar  coordinates  write  x  G  Hn  as  x  —  (u>,r),  where  r  =  d(x,o)  and  oj  G 
Sn_1.  The  hyperbolic  metric  is  then  expressed  by 


ds2  =  dr 2  +  sinh2  r  da;2, 


where  du>2  is  the  usual  metric  in  Sn  1.  Correspondingly,  the  ( n  —  l)-dimensional  area  of 
a  geodesic  sphere  of  radius  r  is 


(2.1)  A „(r)  =  finsinhn  1  r,  where  f2n  = 

The  Laplace-Beltrami  operator  on  Hn  is 

(1  ~  INI2)"  3 


27Tn/2 

rW2) 


is  the  Euclidean  area  of  S 


n—1 


Ah  = 


dx  '  i 
j= 1  L 


(1-Ik||2) 


2\2  — n  ^ 


dx  j  J 5 


which  specializes  to  A u  =  (1  —  \z\2)2d2 /dzdz  (with  z  —  x i  +  ix 2)  in  the  case  of  n  =  2.  In 
polar  coordinates 

3^  3 

A H  =  7-7  +  (n  -  1)  cothr—  +  sinh-2  r  As, 

Or 1  or 

where  As  is  the  Laplace-Beltrami  operator  on  Sn-1. 
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The  space  P(Hn)  denotes  as  usual  the  space  of  all  C°°  functions  with  compact  support 
in  Hn,  i.e.,  it  coincides  with  the  space  X>(B")  of  C°°  functions  /  on  Rn  whose  support 
supp  /  is  in  Bn.  The  Schwartz  space  5(Hn)  of  fast  decreasing  functions  in  H"  is  the  space 
of  C°°  functions  /  on  H"  such  that  for  any  positive  integers  m,  k  we  have 

sup  \AkHf{x)\emd{x’o)  <  oo 
xeH" 


(cf.  [H2]).  This  is  equivalent  to  the  condition  that  for  every  multiindex  a  €  Nn  (where  N 
is  the  set  of  nonnegative  integers),  the  function  Bn  3xh  d^f  /  dxa{x)  has  a  continuous 
extension  to  the  closed  ball  Bn  which  vanishes  of  infinite  order  in  <9Bn.  In  other  words, 
the  space  5(Hn)  coincides  with  the  space  P(Bn)  of  C°°  functions  /  in  Rn  such  that 
supp  /  C  Bn.  The  two  spaces  coincide  even  topologically. 

Fix  k,  with  1  <  k  <  n:  the  space  T  =  Tn  ^  of  totally  geodesic  ^-dimensional  submanifolds 
(k- geodesies  for  short)  of  Hn  is  a  homogeneous  space  under  an  action  of  the  group  50(1,  n) 
of  isometries  of  Hn.  Each  7  6  T  carries  the  ^'-dimensional  area  element  dm //7  induced 
by  the  volume  element  dm  a  in  Hn.  Hence  the  totally  geodesic  ^'-dimensional  Radon 
transform  Rh  —  RH,n,k  is  defined  on  the  space  5(Hn)  by 

RhK 7)-  [  f(x)dmHn(x)  for  all  7  6  F. 

J 1 

The  family  Fx  of  elements  of  V  passing  through  a  fixed  point  x  is  a  homogeneous  space 
for  the  isotropy  group  50(1, n)x  of  x,  which  is  isomorphic  to  50(n).  Hence  carries  a 
normalized  measure  dmriX  which  is  invariant  under  50(1,  n)x,  and  is  ‘independent’  of  x 
in  an  obvious  sense.  For  a  continuous  function  on  T  we  can  define  the  backprojection 
operator  R*H  by 

R*H(f)(x)  =  f  ^(7)  dmr,i(7)  =  f  <f)(gh-^)dh  for  all  x  €  H", 

J  r»  ’  Jso(rx) 


where  g  is  a  fixed  element  of  50(1,  n)  such  that  g  •  o  =  x,  while  h  runs  in  50(n)  and 
dh  is  the  normalized  invariant  measure  in  50(n).  One  of  the  uses  of  the  backprojection 
operator  is  to  find  an  inversion  formula  for  the  Radon  transform.  This  is  based  on  the 
fact  that,  denoting  by  dd  the  area  element  on  the  geodesic  sphere  5(x,r)  of  center  x  and 
radius  r,  we  have 


(2.2)  R*HRHf(x)  =  /  f(y)Tl(d(x,y))dmH(y)  =  f  R(r)  f 

J  H"  Jo  [.J  S 


S(x,r) 


f{y)d0(y) 


dr , 


for  a  function  7Z  —  7Zn,k  on  [0, +00)  (cf.  [H2,  Theorem  1.4.5]):  interpreting  R  as  a  radial 
function  on  Hn  through  TZ(x)  =  7£(d(x,o))  (the  same  abuse  of  notation  will  be  extended 
to  all  radial  functions,  with  no  possible  confusion),  we  write  this  integral  as  7Z*  f(x) — note 
that  the  inner  integral  is  not  normalized—:  in  fact  both  %  and  /  can  be  pulled  back  as 
functions  on  the  group  50(1,  n),  convolved  there,  and  the  result,  pushed  to  Hn  again, 
coincides  with  the  middle  term  of  (2.2). 
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The  function  7Z  turns  out  to  be  [Hi] 

K(r)  =  !r  rp)  r' 

In  [BC]  it  is  shown  that  if  S  is  the  operator  of  convolution  associated  with  the  radial 
function 

<S(r)  =  sinh^-"  r  coshr, 

then  for  an  explicit  polynomial  p  of  degree  k  one  has 

p(&)SR*hRh  =  I 

(where  I  is  the  identity  operator),  which  is  the  filtered  backprojection  formula  mentioned  in 
the  introduction.  For  instance  if  ( n,k )  =  (2, 1)  this  formula  reduces  to  —  (47r)-1  SR*hRe  = 
I :  the  kernel  S  can  be  replaced  by  coth  r  — 1  to  obtain  an  integrable  kernel.  Other  inversion 
formulas,  which  do  not  factor  through  R*hRh,  can  be  found  in  [H3]. 

The  ^-dimensional  hyperbolic  Radon  transform  was  in  fact  modeled  upon  the  Euclidean 
one,  which  is  similarly  defined  in  the  space  Rn.  Namely,  if  /  G  <S(Rn)  and  <j>  G  <S(G), 
where  G  —  Gnfc  is  the  Grassmannian  space  of  all  k-planes — viz.,  of  fc-dimensional  affine 
subspaces  of  Rn — ,  then 


ReH.*)  =  /  f{y)dmE,i r(y)  for  all  7T  €  G, 

J  7 r 

R*E<f>(x)  =  /  (f>(n)  dmG,y(n)  =  /  <f>(gh  • -re)  dh  for  all  y  G  Rn, 

JGy  ’  JSO(n) 

where  dmE  n  is  the  fc-dimensional  Lebesgue  measure  on  it,  while  dmc,y  is  the  normalized 
SO(n)-invariant  measure  on  Gy  =  {it  €  G :  y  €  tt}. 

In  order  to  characterize  the  range  of  the  transform  Rh  in  the  hyperbolic  space  Hn 
we  shall  show  that  it  can  be  intertwined  with  the  Euclidean  Radon  transform.  For  that 
purpose  we  need  to  set  up  some  additional  notation.  Let  E  =  E„  *  be  the  set  of  (n  —  k  — 
1  )-circles,  intersections  of  (n  —  fc)-dimensional  subspaces  of  Rn  with  the  unit  sphere  Sn_1. 
Denote  by  Gq  the  space  of  k- planes  that  intersect  B".  Naturally  the  closure  GE  of  Gq  is 
the  compact  subset  of  G  consisting  of  A:-planes  intersecting  Bn. 

One  can  parametrize  G  and  T  in  a  similar  way.  First  take  polar  coordinates  in  Rn,  so  that 
a  point  x  is  given  by  a  pair  (u,  s )  G  Sn_1  x  [0,  oo).  Then  identify  G  with  the  set  of  triples 
(£,  u>,s)  e  3  x  Sn_1  x  [0,  oo)  such  that  u  G  £:  such  a  triple  represents  the  fc-plane  which  is 
orthogonal  at  the  point  (u>,s)  to  the  (n  —  fc)-plane  through  o  generated  by  £  (observe  that 
(tu, s)  is  the  closest  point  of  the  fc-plane  to  o).  In  particular  G#  =  {(£,u?,s)  G  G:  s  <  1}. 
In  a  completely  analogous  way  T  can  be  parametrized  by  the  same  set  of  triples,  and  the 
(n  —  fc)-plane  generated  by  £  in  R"  projects  to  an  ( n  —  fc)-geodesic  through  o  via  the 
exponential  map  at  o.  Note  that  in  both  cases  the  £  coordinate  is  redundant  if  k  =  n  —  1 
(it  must  be  £  =  {±tu}). 

In  next  lemma  we  find  the  explicit  expressions  in  such  coordinates  of  the  measures  on 
R",  H",  G,  T.  We  use  the  following  notation:  if  tt  G  G  let  dnu>  be  the  fc-dimensional 
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measure  on  {c o  G  Sn_1 :  (w,s)  G  7r};  let  ef(£,w')  denote  the  normalized  rotation-invariant 
measure  on  the  manifold  P  ~  {(£,w')  £  Hx  Sn_1 :  to'  G  £}  (notice  that  dimP  =  (n  — 
k)(k  +  1)  —  1);  finally,  if  y  G  Rn,  let  dy(£,u’)  be  the  [(n  —  fc)fc]-dimensional  measure  on 
{(£,u/)  G  P:  (£,w',s')  G  Gy}  (recall  that  dim  Gy  =  dimH  =  (n  —  k)k).  Similarly  define 
d7w  and  dx(^,u')  in  the  hyperbolic  case. 

LEMMA  2.1.  In  the  coordinates  introduced  above,  the  measures  in  the  spaces  in  consider¬ 
ation  are 


dmg(u>,  s)  =  sn  1  ds  du> 
dm//(w,  r)  =  sinh”-1  r  dr  du 
dmG(£,u',  s')  =  (s,)n_fc_1  ds'  d(£,u') 
dmr(£,uj' ,r')  =  cosh k  r'  sinh n~k~1  r'  dr'  d(^,u') 


(in  Rn), 
(in  Hn), 
(in  G), 
(in  T). 


Moreover  we  have 


dm£iK(u,  s)  =  (sk+1  /s')d„u  (on  tv  =  (^,  u o', s')  G  G), 

sinhfc+1 7* 

dmH^(u,r)  =  v-—-  d7w  (on  7  =  (£,c o' ,r')  G  T), 

dmGiy(Z,u',s')  =  A,(£,w')  (in  Gy  for  y  ~  (w,s)  G  RnJ, 

cosh*"*’1  r; 

dmr,x(^w',r')  =  coshr  <**(£»  ^0  (in  for  x  =  (w,r)  G  Hn). 

PROOF:  The  expression  of  dm#  is  well-known,  while  that  of  drriH  follows  from  (2.1). 

Let  y'  =  (u>',s')  and  x'  =  (w',r'),  the  closest  points  of  7r,  respectively  7,  to  o.  If  /?  is 
the  angle  y'yo  between  straight  line  segments,  then  sin/?  =  s' /s;  therefore  dms,n  equals 
sk /  sin/?  =  sk+1/s'.  Analogously,  if  a  is  the  angle  x'xo  between  geodesic  line  segments, 
then  by  hyperbolic  trigonometry  (cf.,  e.g.,  [C])  we  have  sin  a  =  sinh  r'/  sinhr,  and  obtain 
dm  H,~f. 

To  proceed  with  the  remaining  measures  in  H"  we  must  make  a  preliminary  computa¬ 
tion.  Let  Q  be  the  angle  between  the  geodesic  segment  x'o  and  the  normal  to  7  at  x,  and  let 
A w  =  w'w  be  the  angle  between  w  and  w'.  Since  9  is  complementary  of  the  angle  a,  by  hy¬ 
perbolic  trigonometry  we  have  sin  9/  sin  Aw  =  coshr'  and  cos  9/  cosAw  =  coshr'/  coshr, 
therefore 

d9  cosh2  r' 

^  dK>  =  coshr  ’ 

this  derivative  relates  the  variation  of  the  normal  to  7  at  a  fixed  x — recall  that  the  measure 
dmT  x  is  invariant  under  50(1,  n)x  =  SO{n ) — with  that  of  cP.  (In  the  Euclidean  case  the 
angle  9  simply  equals  Aw,  so  that  the  quotient  in  (2.3)  just  evaluates  as  1.) 

The  factor  (s')n_fc_1  in  dmG  corresponds  to  the  variations  in  the  w'  coordinate  of  tt  = 
(£,w',s'):  the  exponent  is  the  dimension  of  £.  The  same  argument  justifies  the  factor 
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sinhn_<'~1  r'  in  dmr-  To  account  for  the  disparity  between  Au;  and  6  we  let  r1  =  r  in  (2.3) 
and  obtain  coshr',  which  raised  to  the  dimension  of  7  gives  the  other  factor  in  dmp:  this 
corresponds  to  the  variations  of  7  within  the  (k  +  l)-geodesic  containing  7  and  o,  keeping 
the  point  x  fixed. 

The  measure  dmGj!/  equals  dy(£,u')  by  translation  invariance.  The  factor  in  dmp  x  is  the 
product  of  cosn  r'/  coshr,  given  by  (2.3),  and  dim 7  —  1  factors  coshr',  which  correspond 
to  the  variation  of  the  normal  to  7  at  x  in  the  orthogonal  space  to  the  2-geodesic  containing 
o,  x',  x.  | 

With  the  above  described  measures  the  dual  Radon  transform  R*  is  (in  both  the  Eu¬ 
clidean  and  the  hyperbolic  case)  indeed  the  adjoint  of  the  transform  R:  if  /  6  V(Itn)  and 
<f>  £  'D(G)  we  have 


(2.4) 


(ReJ,*) 


f(y)dmE<ir(y ) 


<K*)  dmG(rc)  =  /  f(y)  dmv(y,  tt) 

Jv 


’  Rn 


f(y) 


/  0(tt)  dmG>y(7 r)  dmR(y)  =  (/,  R*E(f>), 

Ugv 


where  dmv(y,n)  =  ds  [sn_1  du>]  dy((,ui')  =  [(s')n~fc~1  d(^,u;')]  [sk  dnuj]  (s/s1)  ds'  is  a  mea¬ 
sure  on  the  manifold  V  =  {(y,7r)  €  Rn  x  G:  y  €  7r)  =  {(o»,s;  s')  G  Sn_1  x  (0,oo)  x 

P  x  [0,oo):  s'  =  (u',u>)s,  Span(o;',a!)  1  £}  (here  Span(u?',a;)  denotes  the  linear  span  of 
a)',  u):  by  [H2,  §1.3.2],  the  measures  in  Lemma  2.1  are  characterized  by  (2.4).  A  similar 
computation  can  be  carried  out  for  Hn. 

3.  Characterizations  of  the  range 

In  the  previous  section  we  saw  that  the  space  S(Hn)  could  be  identified  to  X>(Bn),  and 
D(Hn)  to  £>(13").  In  either  case,  a  function  /  in  S(Hn)  (respectively  Z>(Hn))  can  be 
identified  to  a  C°°  function  in  Rn  with  supp  /  C  Bn  (respectively  supp  /  C  Bn),  and  it 
makes  sense  to  consider  also  the  effect  on  /  of  the  Euclidean  Radon  transform  Re,  in  which 
case  only  the  integrals  along  elements  of  Gg  can  fail  vanishing.  We  then  have  (see  [H2]) 
the  maps 


Rh : 
Re- 


P(Hn)  ->  D(r), 

5(H")  -  S(T), 

£*(Bn)  ->  V(Gb)> 

V(Bn)  ^  V(Gb)  C  V{G). 


Before  we  describe  the  intertwining  operators  between  Rh  and  RE  we  develop  some 
auxiliary  geometric  results.  Let  rj:  Hn  — >  Hn  be  the  homothety  of  ratio  2  with  respect  to 
o,  that  is, 

ry(u;,r)  =  (w,2r)  for  all  (uqr)  £  Hn, 

in  particular 

d(o,r}(x))  =  2d(o,  x); 
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the  map  rj  corresponds,  via  the  exponential  map,  to  multiplication  by  a  factor  2  in  the 
tangent  space  of  Hn  at  o.  If  we  let  i:  Hn  — >  Bn  be  the  identity  map  (which  changes  the 
metric  from  hyperbolic  to  Euclidean),  so  that 

i(w,r)  =  (u;,  tanh(r/2))  for  all  (u>,r)  £  Hn, 
the  map  77  can  be  read  in  Euclidean  polar  coordinates  as 

torjo  i_1(o;,s)  =  (w,tanh(2arctanhs))  =  (w,2s/(l  +  s2))  for  all  (w,a)  €  Bn. 


The  map 


r  =  7077:  Hn  ->  Bn 


is  a  diffeomorphism  that  satisfies 


r(u>,r)  =  (cj,  tanhr)  for  all  (u>,r)  6  Hn; 


we  shall  see  that  r  induces  an  isometry  between  the  conformal  (the  one  we  are  using)  and 
Beltrami’s  model  of  hyperbolic  space.  This  is  shown  by  the  following  proposition: 

PROPOSITION  3.1.  The  image  through  t  of  a  k-geodesic  of  H"  is  the  intersection  with  B" 
of  a  k-plane  of  RT 

PROOF:  Since  r  preserves  the  tangential  component  u  £  S"-1  we  can  assume  n  =  2, 
k  =  1.  Fixed  a  geodesic  7  =  (£',u/,r')  (the  £'  component  is  useless  here,  as  mentioned 
earlier),  its  closest  point  to  o  is  x'  =  (c </,r');  let  x  =  (u>,r)  be  another  point  of  7.  Let  a 
be  the  angle  x'ox  —  and  observe  that  the  angle  xx'o  between  geodesic  segments  is  a 
right  angle. 

In  a  right  hyperbolic  triangle,  the  cosine  of  a  non-right  angle  equals  the  ratio  of  hyper¬ 
bolic  tangents  of  the  lengths  of  the  adjacent  sides  (cf.  [C,  12.99]).  So 

tanhr'  =  tanhr  cos  a. 


Let  y'  —  t(x')  =  (a/, tanhr'):  then  ||y'||  =  tanhr'.  Similarly  y  =  r(r)  =  (w, tanhr)  and 
||y||  =  tanhr.  Moreover  the  angle  y'oy  equals  ui'uj  =  a.  Therefore 

lly'll  =  l!y||cosa 

by  the  previous  formula,  and  from  Euclidean  geometry  one  gathers  that  the  angle  yy'o, 
between  straight  line  segments,  must  be  a  right  angle.  Hence  y  varies  on  the  intersection 
with  Bn  of  the  straight  line  orthogonal  at  y'  to  the  segment  y'o.  | 

REMARK  3.2.  (The  following  description  of  models  for  H"  can  be  found  in  [Y]  ;  see  also 
[Kl].)  Consider  in  Rn+1  the  signature  (l,n)  Lorentz  metric 

y(r',y')  =  x0y0  -  (x,y),  for  x,y  £  Rn  and  x'  =  (x0,x),  y'  =  (yo,y), 
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where  ( • ,  • )  is  the  standard  scalar  product  in  Rn.  The  hyperboloid  model  for  Hn  is  given 
by  the  ‘upper’  sheet 

Q  =  {y(x',x')  —  1  and  x0  >  0} 

of  the  two-sheeted  hyperboloid  {y(x',x')  —  1},  endowed  with  the  distance 
(3.1)  d(x',y')  =  arc  cosh  g(x',y'). 

The  group  of  automorphisms  of  is  thus  50(1,  n),  which  leaves  Q  invariant,  and  the 
k- geodesics  are  the  intersections  with  Q  of  ( k  +  l)-planes  of  Rn+1. 

Beltrami’s  model  for  Hn  is  obtained  by  identifying  Bn  with  D\  =  {x0  =  1  and  ||x||  <  1} 
in  Rn+1,  and  projecting  Q  from  the  origin  of  Rn+1:  a  point  x'  E  Q  gets  mapped  to 
x' /xq  E  D i,  i.e.,  to  x/x0  €  Bn.  It  is  immediate  that  the  fc-geodesics  are  now  intersections 
with  Bn  of  fc-planes  of  Rn.  The  distance  in  D\  is  obtained  by  first  making  the  argument 
of  the  right-hand  side  of  (3.1)  homogeneous,  then  restricting  it  to  Dp.  so 

dix.y )  =  arc  cosh  — (where  x'  =  (1,  x),  y'  =  (1,  y)) 

=  arc  cosh 

As  explained,  e.g.,  in  [C],  or  [M],  identifying  Bn  with  D0  =  {x0  =  0  and  ||x||  <  1}, 
the  composition  k  of  the  inverse  stereographic  projection  of  D0  from  one  pole  of  the  unit 
sphere  {||x'||  =  1}  in  R"+1  onto  the  opposite  hemisphere,  with  the  orthogonal  (‘vertical’) 
projection  onto  D o,  provides  an  isometry  of  the  conformal  model  onto  Beltrami’s.  The  map 
k  is  such  that  the  image  of  each  point  lies  on  the  same  radius  starting  from  o:  therefore 
k  coincides  with  r,  that  is,  the  composition  of  projections  described  above  can  also  be 
obtained  as  a  homothety  in  the  conformal  model  (followed  by  a  change  of  metric).  | 

As  a  consequence  of  Proposition  3.1,  the  map  r:  Hn  — >  Bn  naturally  induces  a  diffeo- 
morphism 

Tk:r-+GB. 

Explicitly 

rfc(£)w ,r)  =  (£,u;,tanhr)  for  all  (£,u>,r)  E  T. 

LEMMA  3.3.  The  jacobian  ofr  along  the  k-geodesic  7  =  (£ ,u>',r ')  at  the  point  x  =  (u;,r)  E 
7  is  coshr'/  coshfc+1  r. 

PROOF:  From  Lemma  2.1,  setting  s  =  tanhr  and  ir  =  77(7),  the  desired  jacobian  is 
obtained  as  the  ratio  of  dm tanhr)  and  dmH,i(w,r).  I 

We  are  now  ready  to  define  the  intertwining  operators  between  Ru  and  Re- 

Proposition  3.4.  For  1  <  k  <  n  let  p:  H"  — >  R  and  o\  T  — ♦  R  be  defined  by: 

p(u,  r)  =  cosh*+1  r,  hence  /»ot'1(w,5)  =  (1-s2)'^+1^2; 
o(£,ut  r)  =  coshr,  hence  a  o  —  (1  —  s2)"1^2. 
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Then  the  transformations  <f> :  D(Hn)  — >  T>(Bn)  and  :  T(r)  — >  T>(Gb)  given  by 

Hf)  =  ( pf)°r~\ 

W)  =  (<*<!>)  OTki 

axe  topological  isomorphisms,  and  the  diagram 


D(Hn) 

Rh 
- > 

®(D 

$ 

$ 

Re 

- ► 

V(Bn) 

V(Gb) 

commutes,  i.e.,  ^Rh  =  Re$- 

(Note  that  p  is  independent  of  n,  and  a  is  independent  of  n,  k.) 

PROOF:  If  7  €  r  and  tt  =  77.(7),  and  if  /  €  X>(Hn),  applying  Lemma  3.3  we  have 

^RhI( tt)  =  <7(7)  f  fdmH, 7  =  Ap/)  0  r-1  =  i?E$/(jr).  I 

Clearly  the  same  proof  furnishes  the  following  result: 

PROPOSITION  3.5.  The  transformations  $  and  \I/  defined  in  Proposition  3.4  also  render 
commutative  the  following  diagram: 


5(Hn)  ■ 

Rh 

- ♦  S{T) 

(3.2) 

$ 

X>(Bn) 

_ 

- *  V(Gb ) 

( and  are  topological  isomorphisms  also  between  these  spaces). 

PROOF:  The  only  observation  that  needs  to  be  made  is  that  $  maps  <S(Hn)  into  V(Bn). 
In  fact  it  is  clear  that  since 

S( Hn)  =  {/  €  C°°(H") :  sup  \Daf(x) |  cosh’71  d(x,o)  <  00  for  every  a  €  Nn,  m  €  N} 

x€H" 

we  have 

^(^(H71))  =  {g  €  C°°(Bn):  sup  |-Day(y)|(l  —  ||y||)-m  <  00  for  every  a  €  Nn,  m  €  N} 

!/€B” 

=  V(Bn).  | 

(The  term  ‘intertwining’  used  here  is  a  slight  abuse  of  language,  since  the  operators  <I> 
and  *3/  do  not  preserve  the  group  actions.) 
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REMARK  3.6.  As  an  immediate  consequence  of  Proposition  3.4  and  Proposition  3.5  we 
obtain  that  any  inversion  of  the  Euclidean  Radon  transform  Re  automatically  yields  a 
corresponding  one  for  i?#,  by  the  relation  RJ{1  —  $~1Rg1' h.  | 

In  a  very  similar  fashion  one  verifies  that  the  map  :  Lp(Hri)  — >  Lp(Bn)  given  by 

$P(f)  =  ( Ppf)°T~\ 

where  pp:  Hn  — >  R  is  the  function 

pp(u,r)  —  cosh^"+1^p  r, 


is  an  isometry  for  every  p  >  1:  in  particular,  the  map  $  itself  is  an  isometry  of  Lp- spaces 
for 


n  +  1 


P  = 


k  +  1' 


Dualizing  the  diagram  (3.2)  we  obtain 


1Rh 


S‘(T) 

- >  S'(H") 

(3.3) 

* 

'Re 

€\Gb) 

- >  £'{ B") 

The  space  <S(T)  can  be  embedded  into  S'(r),  and  1Rh\s(T)  =  whereas  if  restrs  is  the 
restriction  operator  from  Rn  to  B"  then  1Re\d(gb)  =  restr^  oR*E:  however,  by  abuse  of 
notation  the  indication  of  the  map  restr^  will  be  henceforth  suppressed.  We  thus  obtain 
the  diagram 


s(r) 


S'n£(Hn) 


V{Gb) 


£(B") 


In  (3.3)  the  maps  are  isomorphisms.  We  shall  now  proceed  to  make  their  inverses 

explicit  when  acting  on  smooth  functions. 

Lemma  3.7.  The  jacobian  of  r  at  x  —  (w,r)  €  Hn  is  cosh-"-1  r,  and  that  of  Jacobian  of 
Tk  at  7  =  (£,u/,r')  £  T  is  cosh-”-1  r'. 


PROOF:  Same  as  for  Lemma  3.3,  using  Lemma  2.1.  | 

Proposition  3.8.  For  1  <  k  <  n  let  p'\  Hn  -+  R  and  a':  T  ->  R  be  defined  by: 


p'((x>,  r)  —  coshr,  hence  p'  o  r  x(u;,  s)  —  (1  —  s2)  1^2; 
cr'(£,u;,r)  =  cosh*+1  r,  hence  cr'  o  s)  =  (1  —  s2)~^k+1^2. 
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Then  the  transformations  <P':  S'  0  £(Hn)  — >  £(Bn)  and  \P' :  <S(r)  — ►  V(Gb )  given  by 

$'(/)  =  (p7)o'f_1, 

%'{<!>)  =  (<7'(f>)0Tk1 

are  topological  isomorphisms,  and  the  diagram 

R* 

S(T)  - 5'  n  5(Hn) 

*'j 

_  fT£  _ 

V(Gb)  - ♦  £(Bn) 

commutes,  i.e.,  R*H  —  R*^' .  In  other  words,  =  (t(P)-1  and  —  (1'P)~1. 

PROOr:  Recalling  the  value  of  the  jacobian  of  r  from  Lemma  3.7  one  has 

{$/,$'/'}=/  (pfp'f')oT-1  =  f  //'  =  </,/'), 

J  B"  ./H" 

and  since  $  is  an  isomorphism  we  have  (<P*)~1  =  $'  in  the  dense  subspace  S'  H  £(Hn). 
With  VP,  we  proceed  analogously,  and  conclude  that  \P'  =  (’P*)-1,  which  proves  the 
desired  commutativity.  I 

Observe  that  <P,  \P  are  never  unitary:  in  fact  they  do  not  coincide  with  \P',  respec¬ 
tively,  for  any  n,  k. 

For  functions  g  €  X>(R")  Helgason  proved  a  support  theorem  [H2,  Corollary  1.2.25], 
which  characterizes  the  radius  s0  of  the  smallest  ball  in  Rn  centered  at  o  that  contains  the 
support  of  g :  all  is  needed  is  that  J?£<; 1  vanish  at  every  fc-plane  ir  —  (< £,u,s)  with  distance 
s  to  the  origin  greater  than  s0.  Therefore  one  can  see  that  among  the  functions  ip  €  S(G) 
that  are  in  the  range  of  Re  one  can  easily  distinguish  those  that  are  in  RE(T>(3n))  or  in 
R£('D(Bn)),  by  the  conditions  (respectively) 

supp ip  C  {(£,tu,s)  €  G:  s  <  1}, 
supp  ip  C  {(£,w,s)  €  G:  s  <  1). 

Simply  notice  that,  under  either  of  these  support  conditions,  automatically  ip  €  V(G).  As 
to  functions  on  Hn,  the  support  theorem  is  also  valid  for  functions  in  <S(Hn)  (cf.  [H2, 
Theorem  1.4.2]). 

One  way  to  characterize  the  range  of  the  Euclidean  Radon  transform  for  5(Rn)  or 
"P(Rn),  in  case  1  <  k  <  n  —  1,  has  been  to  do  it  in  terms  of  the  solutions  of  a  homogeneotis 
system  of  linear  differential  equations.  The  idea  of  using  differential  equations  originates 
with  John  [J],  who  found  that  the  ultrahyperbolic  operator  in  R4  could  be  used  to  charac¬ 
terize  the  range  of  the  X-ray  transform  (i.e.,  k  =  1)  in  R3.  In  general,  Richter  [R]  exhibited 
a  system  of  translation  invariant,  globally  defined  second-order  differential  operators  which 
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characterize  Rg(S(Tkn))  for  each  k  <  n  —  1.  Recently,  other  characterizations  were  pro¬ 
vided  by  Kurusa  [Kul]  with  a  system  of  ultrahyperbolic  equations,  and  by  Gonzalez  [Gol] 
with  a  single  fourth-order  differential  operator,  invariant  by  all  isometries  of  R".  We  now 
derive  the  corresponding  results  in  <S(Hn)  (and  P(H")).  We  will  concentrate  on  Richter’s 
and  Gonzalez’  operators  and  explain  how  they  give  rise  to  (second-,  respectively)  fourth- 
order  equations  characterizing  the  ranges  of  the  hyperbolic  Radon  transforms,  although 
the  same  procedure  can  be  applied  to  other  differential  equations,  such  as  Kurusa’s. 

Let  E  —  En  be  the  Lie  group  of  isometries  of  Rn,  let  e  be  its  Lie  algebra,  and  let  A 
be  the  universal  enveloping  algebra  of  e.  Denote  by  0  the  left  action  of  E  on  G.  We  say 
that  there  is  a  germ  (at  the  identity)  of  left  action  0  of  E  on  the  open  set  Gb  of  G,  by 
restriction:  for  every  n  €  Gb,  only  those  b  £  E  are  allowed  such  that  0(b) n  £  Gg — these 
form  an  open  set  of  E  that  depends  on  7r  and  contains  the  identity.  Denote  by  v  the  left 
regular  representation  of  E  on  <S(G),  so  that  v(b)ip  =  ip  o  0(b)~l  for  every  if  £  S(G ),  and 
extend  to  A  the  infinitesimal  left  regular  representation  dv  on  S(G )  in  the  usual  way:  thus 

dm 

du(V\  •  •  •  Vm)ip(ir)  =  — - x— {^(expfiVi  •••exptmym)^(7r)]fl=...=:tm=o 

at\-  ■■  at m 

dm 

=  dt\  T."  dt . W(exP(-Dnbm)  •  •  •  exp(-^y1))7r)]tl=...=tm=o 

for  each  Vi, . . . ,  Vm  €  e,  ip  £  «f>(G),  it  £  G. 

The  expression  after  the  last  equality  sign  can  be  used  to  define  the  ‘infinitesimal  left 
regular  representation’  dv  on  V(Gb)  (even  though  the  ‘finite’  left  regulax  representation 
v  on  V(Gb)  does  not  itself  make  sense).  In  fact,  on  each  fixed  n  £  Gg  the  element 
exp(— t mVm)  •  •  ■  exp(— t\Vi)  of  E  can  act  whenever  —a  <  t\, . . .  ,tm  <  a,  for  some  positive 
a  depending  on  7r,  Vi,  . . .  ,Vm. 

Conjugating  0  by  77,  we  get  a  germ  of  left  action  a  of  E  on  T:  that  is,  a(b)  ==  rf1  o 
0(b)  o  T)..  Notice  that  a(b )  is  an  automorphism  of  Hn  if  and  only  if  it  keeps  o  fixed.  (The 
other  elements  of  E  cannot  act,  through  a,  as  hyperbolic  isometries,  since  they  Eire  only 
defined  on  part  of  T,  and  push  the  remaining  ^-geodesics  ‘off’  Hn.)  By  analogy  to  the 
above  situation,  define  the  infinitesimal  left  regulax  representation  dfj.  on  <5 (IIn)  by 

dm 

dn(Vi  ■  •  •  Vm)<f>( 7)  =  -  [6(a(exp(-tmVm)  ■  ■  •  exp(-fi  Vi))7)]tl=...=*m=0 

for  each  Vx, . . . ,  Vm  £  e,  f>  £  S(T),  7  €  T. 


Of  course  dfx(V\  ■  ■  ■  Vm)<p  =  [dv(Vi  •  ■  ■  Vm)((j>  oq'jjo  Tfc. 

Let  Tj  £  e  be  the  ‘infinitesimal  translation  in  the  j- th  coordinate’  of  R",  for  j  —  1, . . . ,  n: 

identifying  e  with  the  Lie  algebra  of  matrices  f  ^  ^  V  where  L  £  so(n)  and  v  £  Rn,  then 


Tj  corresponds  to  the  matrix  all  of  whose  entries  all  vanish  except  the  (j,n  +  1),  which 
equals  1.  Also,  let  X; j  be  the  ‘infinitesimal  rotation  around  the  origin  in  the  plane  of 
the  i-th  and  j-th  coordinates’,  for  qj  =  l,...,n  with  i  ^  j:  the  only  nonzero  entries 
of  its  matrix  are  the  (i,j),  which  equals  1,  and  the  (j,i),  which  equals  — 1.  Richter’s 
‘pre-operators’  V,ji  £  A,  for  distinct  i,j,  k  =  1, . . .  ,n,  are  given  by  [R] 


Viji  —  TiXji  -f  TjXii  +  TiXt] , 


12 


whereas  Gonzalez’  is  [Gol] 


v'=  E  w 

l<j<j</<n 

Richter’s  differential  operators  du(Viji)  on  S(G)  are  invariant  by  translations,  but  not 
by  rotations,  while  Gonzalez’  operator  dv(V)  enjoys  both  invariances.  Notice  that  Viji 
and  V  are  independent  of  k,  whereas  the  differential  operators  proper  are  not,  as  they  act 
on  functions  defined  on  different  Grassmannian  manifolds.  Stipulating  that  a  differential 
operator  D  on  S(Gb)  is  said  to  be  invariant  by  a  motion  b  €  E,  such  that  0(b)G b  intersects 
Gb,  if 

D[tp  o  0(b)]  =  (Dip)  o  0(b)  on  {x  €  Gq  :  0(b) rt  6  G b)  for  every  ip  6  S(G b) 

(the  equality  would  mean  that  D  commutes  with  v(b),  if  the  latter  made  sense),  the  second- 
order  operators  D,ji  =  du(Viji)  on  S(Gb)  are,  as  well,  invariant  by  translations  but  not 
by  rotations,  while  the  fourth-order  D  —  du(V )  is  fully  invariant  under  E. 

Pulling  back  with  T,  it  is  easy  to  see  that  the  differential  operators  Ciji  —  on 

S(r)  are  given  by 

(3.4)  Cijip  —  a~ldfi(Viji)(o<p)  for  every  p  6  «S(P).  ; 

It  is  easy  to  realize  that  CVp  are  neither  invariant  by  rotations  around  o  (becairse  D{ji  are 
not),  nor  by  other  automorphisms  of  Hn;  yet,  they  are  invariant  under  the  germ  of  action 
a  of  the  subgroup  of  E  of  Euclidean  translations  in  Rn. 

Analogously  pulls  back  D  to  the  fourth-order  operator  C  —  on  <S(P),  given 

by  an  expression  similar  to  (3.4).  Now,  C  is  fully  invariant  under  the  germ  of  action  a  of 
E ,  whereas  the  only  hyperbolic  isometries  which  leave  it  invariant  are  those  obtainable  as 
a  of  some  element  of  E,  i.e.,  the  rotations  around  the  origin. 

The  explicit  computation  of  Cpj  and  C  is  relegated  to  the  appendix  to  this  paper.  We 
summarize  the  above  remarks  in  the  following  statement. 

PROPOSITION  3.9.  If  1  <  k  <  n  —  1,  a  function  <p  in  N(r)  (or  D(T))  is  in  the  range 
RhS( Hn)  (respectively  RHV(Hn))  of  the  k-dimensional  hyperbolic  Radon  transform  if 
and  only  if  it  is  annihilated  by  the  second-order  differential  operator  Ciji  for  all  1  <  i  < 
j  <  l  <  n;  equivalently,  if  and  only  if  it  is  annihilated  by  the  fourth  order  differential 
operator  C .  | 

The  remaining  case,  namely  k  ~  n  -  1,  is  covered  by  another  kind  of  characterization. 
Namely,  for  all  1  <  k  <  n  —  1,  the  following  Euclidean  moment  conditions  characterize  the 
functions  ip  €  D(G)  which  are  in  the  range  of  the  Radon  transform  Re  acting  on  X>(Rn) 
[H2,  Corollary  1.2.28]: 

For  every  rn  6  N  there  exists  a  homogeneous  polynomial  Pm  on  R"  of  degree  m  such  that 
f  f  sm+n  k  1  ip(^,u,s)ds  du>  —  Pm(ui')  for  each  £  £  E  and  u/  € 

v  ,v  0 

(Such  conditions  also  characterize  the  range  for  functions  of  5(Rn),  but  only  in  the  case 
k  —  n  —  1.) 

The  corresponding  result  in  the  hyperbolic  space  is  the  following: 
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THEOREM  3.10.  Let  <j>  £  L>(T),  respectively  «S(T).  Then  necessary  and  sufficient  condition 
for  <j>  to  belong  to  R^{V{ Hn)),  respectively  R^(S( Hn)),  is  that  for  every  m  £  N  there 
exists  a  homogeneous  polynomial  Pm  on  Rn  of  degree  m  such  that 


tanh 


m-fn-  k  —  1 


coshr 


<f>(£,u,r)dr 


du  —  Pm(u') 

for  each  £  £  E  and  u>'  £  £. 


These  are  called  the  hyperbolic  moment  conditions. 

PROOF:  Recall  that  T(<^)  belongs  to  D(B"),  respectively  P(Bn).  Hence  a  necessary  and 
sufficient  condition  for  ^(4>)  to  be  in  the  range  of  the  Euclidean  Radon  transform  R j?  is 
that  T(</>)  satisfy  the  Euclidean  moment  conditions.  Substituting  (£,a>,s)  by  Tk(£,u),r)  in 
them  one  obtains  a  g  in  P(Bn),  respectively  P(Bn),  such  that  =  Re{q)  if  and  only 
if  <f>  satisfies  the  hyperbolic  moment  conditions.  Let  /  =  then  /  is  in  D(Hn), 

respectively  <S(Hn),  and 

RhU)  =  *-1Re$U)  =  ^~lRE{g)  =  0.  I 


Appendix:  Computation  of  the  range-characterizing 

DIFFERENTIAL  OPERATORS 

Euclidean  case.  For  distinct  indices  a,  b,c  and  distinct  i,j,  l  let  (  —  C(a,b,c,i,j,i ) •'  R4  — >  E 
be  the  map  given  by 

C(*i>*2,<3,<4)  =  exp(-tAXab)exp(-t3Tc)exp(-t2Xij)exp(-t1Ti) 

for  ti,t2,i2,tA  €  R.  Let  {ei,...,en}  be  the  standard  basis  of  Rn.  Let  t  r{j(t )  be  the 
one-parameter  subgroup  in  E  generated  by  —Xij:  so  =  exp (—tXjj)  =  is 

the  rotation  of  angle  t  in  the  (i,j)  plane.  Such  subgroup  has  natural  actions  on  E  and  on 
Sn_1,  which  will  merely  be  denoted  by  juxtaposition.  Let  P$  be  the  orthogonal  projection 
of  Rn  onto  (  6  E,  For  every  fc-plane  7 r  =  (£,w,s)  £  G  and  for  small  in  absolute 

value  we  have 


=  ( rab(ti)rij(t2)^ , 

rab(t4)rij(t2 )  dir(su;  -  -  hP^r ji{t2)ec), 

||sw  -  tiP^i  -  t2P^ji{t2)ec i|), 

where  dir(u)  =  v/  ||u||  for  0  ^  v  £  Rn.  Set 


Ci(tt) 


dt  i 


[C(*i,0,0,0)(7r)]tl=o, 


C234(7r)  — 
Cl234(7r)  = 


d 3 


dt2  dt3  dtA 

cP 


dt i  dt2  dt2  dtA 


[C(0,  ^2)L3)^4)(7r)]<2  =  (3  =  <4=0) 


l((t\,t2,t2,  t4)(7r)]t1=<2=<3=<4=o 
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for  all  possible  combinations  of  order  one  through  four. 

Therefore 

di 

du(TlX{jTcX  ~q^  ~q^  [ip  °  C(^l  5  ^2>  ^3 1  ^4  )] <i  =<2  =  <3  =  <4  =0 

=  d4ip(  Ci ;  C2 ;  (3;  (4) 

+  <^3V,[(Ci;  C2;  C34)  +  (Ci;  Cs;  C24)  +  (Ci;  (4; (23) 

+  (C2;  Ca ;  C14)  +  (C2;  C4 ;  C13)  +  (C3;  (45(12)] 

+  d20[((i2;  (34)  +  ((13;  (24)  +  (Cm;  (23) 

+  ((1 ;  (234)  +  ((2;  (134)  +  ((3;  (124)  +  ((4;  (123)] 

+  dy{C  1234)) 

where,  denoting  by  P  —  P^riu1-  the  orthogonal  projection  of  Rn  onto  (  fl  m2-,  and  by  v  ■  w 
the  ordinary  scalar  product  of  v,w  €  Rn, 

CiO)  =  (0,  - Pei/s ,  -lo  ■  e/), 

(2  (  tt)  =  (-Aby^-ApugO), 

(3(7t)  =  (0,  — P ec/s,  -w  ■  ec), 

(4(tt)  =  (— A'a&(,  —XabU,  0), 

(12(7r)  =  (0,XijPe</5,0), 

(13(71)  =  (0,  -[(w  •  ei)Pec  +  ( u  ■  ec)Pei  +  (e/  •  Pec)u]/s 2,  e,  •  Pec/s ), 

(14(11)  =  (0,  XabPet/s,0), 

(23(71)  =  (0,  [XijP  -  PXij]ec/s,  -uj  •  A,yec), 

(24(71)  =  (XabXij(,XabXijLO,  0), 

(34(71)  =  (0,  XabPec/s,0), 

(123(71)  =  (0,  [(w  •  ec)XijPet  +  (w  •  ei)[XijP  -  PXij]ec  -  ( u  •  Ayec)Pe; 

(e;  ■  P A ,yec)tu  +  (^/  ■  P^c)A ]jCu]/s  ,  e;  •  PA\yec/.s), 

(124(71)  =  (0,  -XabXijPei/s,  0), 

(134(71)  =  (0,  Aa6[(u>  •  ei)Pec  +  (w  •  ec)Pet  +  (e/  •  Pec)u}/s2, 0), 

(234(71)  =  (0,  -Xab[X ijP  -  P A,j]ec/ s, 0), 

(1234(71)  =  (0,  -Arafc[(w  •  ec)A ijPe,  +  (a;  •  e,)[ApP  -  PA,y]ec  -  (uj  •  A',yec)Pe, 

-  (e,  •  PXijec)u  +  (e/  •  Pec)Atyu;]/,s2, 0), 

where  we  stipulated  that 

-Xu* 

XabXijt 


A 

d2 

=  dt2  i=«2=o- 
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We  obtain  the  operator  D  as 

D  =  du(V)  =  Y,  +  dvWiXjiTjXn)  +  du(TtX pT,^) 

i<i<i<i<n  +  du(TjXi,TiXji)  +  dviTjXnTjXu)  +  du^XnTiXa) 

+  duiTtXijTiXj,)  +  dviTtXijTjXu)  +  dv{T\XijTiXij)}. 

As  to  D{ji  we  simply  have 

Dm  =  du(Vijt )  =  du(T,Xjt)  +  dviTjXn)  +  du(T,Xtj) 
for  i,j,  l  =  1, . . . ,  n,  where 

d2 

dv(TiXXJ)i>  =  0  C(*i,<2,0,0)]tl=<2=0  =  d20(C i;  C2)  +  #(C  12), 

and  Ci,  C2,  C12  are  given  above. 

Pullback  to  the  hyperbolic  case.  As  described  in  Section  3,  the  pullback  of  the  oper¬ 
ator  D  is  C  =  so  its  effect  on  4>  e  <S(r)  will  be 

C</>  =  o~lD([o4>\  0  r^1)  0  Tk. 


The  summands  are 

v-'duWXijTcXabW  =  d'^-e  2-ez-ex) 

+  d3</>[(0i;  02;  ^34)  +  (#i!  03;  $24)  +  ($1;  $4;  #23) 

+  ($2!  $3!  $u)  +  {6 2;  $4;  ^13)  +  (^3;  $4i  ^12 )J 

+  d2<£[(0i2;  ^34)  +  (^13^24)  +  (^14^23) 

+  (#1;  $234)  +  ($2!  ^134)  +  ($3!  ^124)  +  (&t;  #123)] 

+  ^(^1234) 

+  tanhr[cfV[(02;0l;04)0[  +  (0l;  #2;  04)03] 

+  d26[(62 ;  ^4 )0\  +  (0£;  e2<)0[  +  (01;  024)03r 
+  W2Z)8\  +  (0(;  04)02r3  +  (02;  014)03r 

+  (02;  04)0i3  +  (04 ;  0i2)03] 

+  ^^[(^24)^13  +  {0u)02  3 

+  (^234)^1  +  (0124)^3  +  (^4)0123]  ] 

+  (1  +  2  tanh2  r)  [d2<j>{62-  04)0[0J 

-bd0[(024)0[03  +(04)0[0^1], 

where  0j,  for  each  multiindex  J  C  {1,2, 3, 4},  is  obtained  from  (j  by  substituting  in  its 
expression  the  variable  s  with  tanhr  and  multiplying  by  cosh2  r  its  third  component,  which 
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is  then  denoted  as  8 and  8'j  is  obtained  from  9 j  by  doubling  its  third  component.  Thus 
for  instance 


6\  —  (0,  —Pei/  tanh  r,  —  u  ■  ei  cosh2  r), 

9\  —  (0,  —Pei/  tanhr,  —2 u  ■  e/  cosh2  r), 

=  —  u>  •  ei  cosh2  r. 

Observe  that  9rj  vanishes  (and  consequently  9'j  —  9 j)  unless 

J  =  {!})  {3},  {1,3},  {2, 3},  {1,2,3}; 

the  above  sum  does  not  appear  symmetric  in  the  indices  1,2, 3, 4  because  the  terms  that 
vanish  were  omitted,  and  because  9 j  has  been  replaced  by  9'j  only  in  case  they  were 
different.  Furthermore  the  factor  Qrj,  when  appearing  in  a  summand,  is  meant  to  multiply 
the  entire  corresponding  differential,  not  its  argument:  thus 

dVKfoi&W +  («!;«*;  <W] 


stands  for 

<J[dV(«2;  M<)  +  ft.)- 

The  whole  operator  C  is  now  obtained  in  the  same  way  as  D  in  the  preceding  section. 
Analogously  we  obtain  the  pullback  Ciji  of  the  operator  D^i  from 

<H-xdu{TiXtJ)^<f)  =  dV(0i;02)  +  <ty(0i2)  +  tanhr0[ch/>(02). 
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